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This paper presents particle swarm optimization (PSO) based approach for 
solving optimal reactive power dispatch for minimizing power losses. The 
control variables are bus voltage magnitudes (continuous type), 
transformer tap settings (discrete type) and reactive power generation of 
capacitor banks (discrete type). The algorithm solution of PSO is tested on 
a standard IEEE 30 Bus system. The intention is to minimize power losses 
by optimizing the reactive power dispatch with optimal setting of control 
variables without violating inequality constraints and satisfying equality 
constraint. The detailed results for different cases have been listed 
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I. INTRODUCTION 

Power system is a very bulk, complex and 
interconnected set-up having generation, transmission, 
distribution and loads. The power system operates in a 
constantly changing environment as generator outputs, 
loads and key operating parameters keep changing 
continually. The focus has been on implementation of 
equipments that can keep the power system reliable. In 
modern time, power networks operate under highly 
stressful situations as power demand is continuously 
increasing so the analysis of power system is don't while 
considering thses factors. The major concerns for power 
systems are: minimization of cost, minimization of losses, 
cost and stability (both voltage stability & power sytem 
stability) [1]. 

Optimization in power system governs with the 
overall economics of the full system. The power system 
optimization consists of best sizing and placement of 
reactive power resources (reactors, capacitors, and SVCs). 
The control settings like active power of generators, 
electrical device faucets, and regular voltages of generators 
are optimized not just for the base case system 
configuration however conjointly for many completely 
different system configurations. 

Load flow or Power flow solution is essential for 
continuous evaluation of power system, planning the 
control and best operation of power system. The quest of 
any power system is to run the system optimally. The 
purpose of an Optimal Power Flow (OPF) algorithm is to 
search out steady stateoperation point that minimize 
generation cost, loss etc. or maximizes financial aid, 
loadability etc. while maintaining satisfactory system 
performance in terms of limits onreal and reactive powers 
of generators, line flow bounds, output of numerous 
compensatingdevices etc. The proper selection and 
coordination of apparatus for governing reactive power 
and voltage stability are amongst the main tasks of power 
system engineers. Conventionally, classical optimization 
methods were used to efficiently solve OPF. In recent 
years, Artificial Intelligence (AI) methods have been 
emerged that can solve extremely advanced OPF issues. 

Reactive power optimization is having significant 
importance for both in day-to-day operations of power 
systems and for future planning. It utilizes all the reactive 
power sources judiciously, while forecasting suitable 
location and size of VAR compensation in a system. The 
financial side of reactive power planning and scheduling 
have a remarkable effect on the profitable and reliable 
operation of a power system as the fuel costs and capital 
investments are increasing day by day [2]. The electric 
power systems all over the world are moving in the 
direction of decontrolled or deregulated electricity 
markets. Additional services like frequency control, 
system control and system restart are required to control 



frequency, security, stability and voltage profile of the 
system and to safeguard the generation and transmission. 
Reactive power and voltage control is a mandatory service 
to sustain voltage profile through injecting or absorbing 
reactive power in electricity market. 

The first comprehensive survey related to optimal 
power dispatch was given by H.H. Happ which reviewed 
the development of optimal dispatch (or economic 
dispatch) and he summarized the different methodologies 
to OPF along with their limitations and with both single 
area and multiple area cases. Afterward an IEEE working 
group presented bibliography survey of major economic- 
security function which was mainly related to operating 
economicsof the system. In the bibliography survey, the 
economic dispatch was conferred under "Real Time 
Operation" function with different categories [3]. 
Subsequently, Carpentier conferred a survey and classified 
the optimal power flow algorithms established on their 
solution methodology which included the power flow 
equations, generator real and reactive power constraints, 
bus voltage magnitude constraints, and bus voltage angle 
difference constraints for buses connected by transmission 
elements. If voltage and angle are taken as variables in 
place of P and Q, the restriction of fixing the reference 
voltage can be lifted.Later, Chowdhary prepared a survey 
report on economic load methods which pointed out the 
importance and related area to economic load dispatch, 
optimal power flow, economic dispatch related to AGC 
etc. [4]. E. Lobato et al. anticipated LP centered OPF for 
optimization (minimization) of transmission losses and 
generator reactive margins of the Spanish power system. 
N. Grudinin suggested a reactive power optimization 
model that was based on successive quadratic 
programming (SQP) methods. A wide spread variety of 
conventional optimization methods have been applied in 
deciphering the RPO problems considering a single 
objective function such as Newton-based techniques [5,6], 
quadratic programming [7], linear programming[8-9], non- 
linear programming [9], Sequential unconstrained 
minimization technique [10], interior point methods [11] 
and parametric method [12]. 

For a power system, the difficulty of reactive power 
planning may be categorized to be an optimization 
problem for which numerous methods have been proposed 
to solve. For solution of all optimization problems, no 
known single optimization method is available. A lot of 
optimization methodologies have been recognized in 
recent years for deciphering different kinds of optimization 
problems. The conventional optimization methodologies 
are: Linear programming (LP), non -linear programming 
and gradient based techniques for solving reactive power 
optimization problems [13-16]. As linearized models use 
approximations, consequently optimal results are not 
signified by LP for objective function being utilized in 
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reactive power optimization problem. Adding to that, 
conventional or Traditional solution strategies have 
tendency to converge to a local optimal solution instead of 
the global one. 

Expert System methodologies have also been 
recommended for reactive power based calculations, which 
are based mostly on 'if-then' dependent rules. The 
evolutionary computational methodologies i.e. 
Evolutionary programming (EP), Genetic algorithm (GA) 
and Evolutionary strategy have moreover been estimated 
to solve the optimizations troubles involving to the 
reactive power dispatch [17-19]. For deciphering complex 
engineering difficulties, the modern (non-traditional) 
optimization approaches are very influential and accepted 
approaches. These approaches are genetic algorithm, 
neural networks, fuzzy optimization, ant colony 
optimization and particle swarm optimization algorithm. 

Particle Swarm Optimizer is a population -based 
stochastic method for global minimization of objective 
functions [20]. Objective functions are a way of 
quantifying everyday real world problems that describe 
properties that need to be minimized to obtain some 
particular outcome. Often objective functions can have 
many parameters that will influence the property that is 
being optimized. Objective functions will have a number 
of characteristics that help determine the how well it will 
be optimized. One of its main attractions is its ability to 
find optimal solutions without the need to compute 
derivatives. It is, in summary, multi-point derivative-free 
optimizer. 

II. Reactive Power Optimization 

Reactive power is essential for reliable operation of 
the bulk power system as it supports power flow. The 
reactive power is an indispensable element of the AC 
transmission grid. The demand of reactive power changes 
at a greater rate than the active power for the same change 
in voltage. Throughout the normal regular operation, 
power systems can go through both over-voltage and 
under-voltage violations which can be overcome by 
voltage or reactive power control. By monitoring the 
production, adsorption and flow of reactive power on all 
the stages in the system, voltage or reactive power control 
can retain the voltage profile inside permissible limit and 
decrease the transmission losses. Generators connected for 
transmission are usually essential for supporting reactive 
power flow. 

Reactive power optimization (RPO) is one of the 
difficult optimization problems in operation and control of 
power system. RPO problem is a special tool to obtain the 
optimal state of the control variables by minimizing the 
certain objectives while satisfying the equality and 
inequality constraints. The most commonly used objective 
is loss minimization. The control or independent variables 



of power system are real power generation excluding lack 
bus, voltage magnitudes at generator buses, transformer 
tap settings and reactive power injection due to 
capacitor/inductor banks. The dependent or state variables 
are load bus voltage magnitudes and angles, slack bus 
power generation, reactive power generation at generator 
buses and transmission line loadings. 

A comparatively modern, new and powerful technique 
Particle Swarm Optimization (PSO) has been practically 
revealed to execute well on numerous optimization 
problems [20-21]. The population based stochastic 
optimization technique i.e. PSO algorithm is applied while 
satisfying equality constraints and not violating inequality 
constraints. The endeavour of minimising reactive power 
losses is achieved by correct adjustment of reactive power 
variables like reactive power generation of capacitor banks 
(Q ci ), generator voltage magnitudes (VgO and transformer 
tap settings (t k ) [12]. 

The recommended PSO algorithm solution is tried on 
the typical IEEE 30-Bus test system with both discrete and 
continuous control variables while keeping the system 
under safe voltage stability limit. The suggested algorithm 
shows improved results. 

III. Particle Swarm Optimization 

Particle Swarm optimization (PSO) is a comparatively 
new evolutionary algorithm that is used to find 
optimum/best (or close to optimal) solutions to qualitative 
and numerical issues. PSO applies the conception of social 
interaction to problem solving. 

Particle swarm optimization was originally established 
by James Kennedy and Russell Eberhart [20]. The 
Eberhart and Kennedy model makes an attempt to seek out 
the most effective compromise between its two main parts, 
individuality and sociality. Particle swarm optimisation 
(PSO), that is a population based stochastic optimization 
technique, shares several similarities with evolutionary 
computation techniques like genetic algorithms (GA). The 
system is initialized with a population of random realistic 
solutions and searches for best by modernizing 
generations. However, unlike GA, PSO has no evolution 
operators like crossover and mutation. PSO algorithmic 
program has additionally been demonstrated to perform 
well on genetic algorithmic program test function. In PSO, 
the potential solutions, referred as particles, fly through the 
problem space by following the present optimum particles. 

In a PSO algorithm, particles modify their positions 
by flying around in multidimensional search area till a 
comparatively unchanged position has been encountered, 
or till machine limitations are exceeded. In scientific 
discipline context, a PSO system combines a social-only 
model and a cognition-only model. A particle changes its 
position using these models. Every particle keeps track of 
its coordinates within the problem space that are related to 
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the most effective solution, fitness; it's achieved to this 
point. The fitness value is additionally stored that is named 
P best . Another best value that is traced by the optimizer is 
the best value, achieved to this point by any particle within 
the neighbors of the particle. This position is named l best . 
Once a particle receipts all the population as its topological 
neighbors, the most effective value could be a global best 
and is named G bes t. At each step PSO changes the velocity 
(accelerates) of every particle toward its P bes t and l best 
locations (local form of PSO). Acceleration is weighted by 
a random term, with discrete random numbers being 
produced for acceleration toward P best and l best locations. In 
past proven years, PSO has been successfully applied in 
research analysis and application areas. 

It's incontestable that PSO gets better results in a 
quicker, cheaper manner compared to alternative methods. 
Another reason that PSO is eye-catching is that there are 
not many parameters to regulate. One version, with slight 
variations, works well in a big range of applications. 
Particle swarm optimization has been utilized for 
approaches which can be used across a large range of 
applications, also as for specific applications centered on a 
specific demand. In the past many years, PSO has been 
successfully applied in several research and application 
areas. 

The fundamental terms employed in PSO technique are 
[23, 24]: 

Particle X (i): It is a candidate solution described by a k- 
dimensional real-valued vector, where k is the no. of 
optimized parameters. 

Population: It is basically a set of n particles at iteration i 
and at time T. 

Swarm: Swarm is defined as an apparently unsystematic 
population of moving particles that tend to bunch together 
while each particle appears to be moving in a random 
direction. 

Particle velocity V(i): Particle velocity is the velocity of 
the moving particles signified by a d-dimensional real 
valued vector. 

Inertia weight w(i): It is a regulation parameter, which is 
used to regulate the impact of the past (previous) velocity 
on the present velocity. Hence, it effects the trade-off 
between the global and local exploration capacities of the 
particles. For the initial stages of the search method, large 
inertia weight to reinforce the global exploration is usually 
recommended while it must be reduced at the last stages 
for higher local exploration. Therefore, the inertia factor 
drops linearly from about 0.9 to 0.4 throughout a run. 
Individual best X*(i): When particles move through the 
search space , it matches its fitness value at the existing 
position to the best fitness value it has ever grasped at any 
iteration up to the current iteration. The best position that 
is related with the best fitness faced so far is called the 
individual best X* (i). 



Global best X (t): Global best is the finest position 
amongst all of the individual best positions achieved so 
far. 

Stopping criteria: Termination of the search process will 
take place whenever optimal control variable settings are 
obtained. 

IV. Problem Formulation 

The objective function 'minimization of real power 
losses', by reactive power optimization, is a constrained 
optimization problem. 

Here IEEE 30 bus system is used for applying PSO 
algorithm. The optimal Reactive Power Dispatch problem 
targets at minimizing the real power loss in a power 
system while satisfying the unit and system constraints. 
The purpose is achieved by correct adjustment of reactive 
power variables like transformer tap settings (T k ), reactive 
power generation of capacitor banks (Q a ) and generator 
voltage magnitudes (V G 0 

The minimization of system real power losses (MW) is 
formulated: 

nl 

F = Min.P^ = ^Skivf + " ^ ™(*«- 'j)] 

E=l (1) 

The real power loss given by P Loss is a non-linear function 
of phase angles and bus voltages that are functions of 
control variables. 

In formulation, 

g k : Conductance of the k th line; 
ni : The number of transmission lines; 
Vi &Vj : The voltage magnitude at the end buses i& j, 
: The voltage phase angles at end buses i& j. 




Fig.l: One Line Diagram of IEEE-30 Bus System 
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The above minimization is subjected to the 
following equality and inequality constraints: - 
Equality constraints: 

These constraints are usual load flow equations that 
are defined as [24]: 

1. Real Power Constraints: 

Pa -Pz>i-Vi feiCtufai ~ Sj) + BySmfa - 

(2) 

i= 1,2 N B 

2. Reactive Power Constraints: 

Qai ~ Qzi ~ P[ ^ Vj ((kjSinfa - 6 } ) - ^-Onfe - 
5,))=0, 



(3) 



i = 1,2 

where, 

Vi 



Gjj.. Bjj 
N B -1 



X 



Real & reactive powers injected into 
network at bus i 

Magnitude of Voltage at bus i 
Magnitude of Voltage at bus j 
Voltage angle difference between bus i and 
busj 

Mutual conductance and susceptance 

between bus i and bus j 

Reactive power generation at bus i 

Total number of buses not including slack 

bus 

Number of PQ buses 



Inequality constraints: 

3. Bus Voltage magnitude constraints: 

v .min ^ v _ ^ y.ma* , [eNq 

4. Transformer Tap position constraints: 

5. Generator bus reactive power constraint: 

6. Reactive power source capacity constraints: 



f L 



7. Transmission line flow constraints: 



8. Generation capacity constraint: 

P^ iD < P &i < ?S SX : i*N B 



(4) 
(5) 
(6) 

(7) 
(8) 

(9) 

(10) 



The total power generation is supposed to cover 
the P D (the total load demand) and P L (the real power loss 



in transmission lines). The relation is regularly expressed 
as: 

9. Power balance constraint: 

^P Gi =Po+PL (ID 

The symbols used are as follows: 
Tap setting of transformer at branch k 
Reactive power generated by E th capacitor bank 
Reactive power generation at bus i 

Apparent power flow through the E th branch 
Number of generator buses 
Number of capacitor banks 
Total number of buses 
Conductance of buses 

Number of tap-setting transformer branches 



t k 

Si 

N g 

N c 

N B 

gk 

X- 



The control variables for voltage- control 
problem, which is able to be modified by the Particle 
Swarm optimization process, are: 

• Voltages magnitude at voltage-controlled buses 
(PV-buses) together with the slack bus. 

• Adjustable shunt capacitor banks. 

• Transformers tap settings. 



V. Methodology 

The various steps concerned with the implementation 
of particle swarm optimization to the reactive power 
optimization problem are [23-24]. 



Step 1: 



Step 2: 



Step 3 \ 

Step 4: 

Step 5: 
Step 6: 



Primarily the input parameters of the system 
(bus, line and generator data) are scanned and 
the lower and upper boundaries of every 
variable are also identified. For N generators, 
optimization is applied for N-l generators and 
generator of maximum capacity is considered at 
slack bus. 

The particles of the population are arbitrarily 
initialized i.e. are arbitrarily selected between 
the corresponding minimum and maximum 
values. Also assign the velocity V firstly 
between [-1 and 1]. 

Achieve power flow solution and calculate 
losses via Newton-Raphson method. 

The best fitness is assigned asP best . At this 
juncture the P best is also the G\, est . 

Iteration /= i+1 is modernized. 

Adjust the inertia weight w given by 



W = 



X iter 



(12) 
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Update the velocity v of each particle according 
to the stated equation. 



Step 7: 



V ft,/, i + D = V ft,;, 0 + C L *rond * (pbsst* ft) - 
X(kJ t zD)-h C 2 Vrand*(£ tof ^(fc) - £ ft,;, 0) 

(13) 

Step 5: Position of each particle is modified too as per 
the stated equation. If a particle disrupts the 
position bounds in any dimension, its position is 
fixed at the right limit. 



*ft,/, i + 1) = xtkj -1,0 + v(kj f 0 



(14) 



Step 9: Assessment of each particle is prepared along 
with its updated position by running power flow 
and evaluate the fitness function. If the 
evaluation/assessment value of every particle is 
improved than the former P best then the present 
value is fixed to be P bes t. If the best P bes t is better 
than G b est, the value is fixed to be G bes t. 

Step 10: If one of the stopping benchmark is fulfilled 
then we go to Step 11. Else, we go to Step 5. 

Step 11: Gbest is the optimal/best value which is newest 
generated by the particle. 



Read the input Data 
Reactive Power (Q), P, V 



I 



Initial Searching Points and 
velocities are randomly generated 
within their limits 

1 



Pbest is set to each initial searching 
point. The best evaluated values 
among Pbest is set to Gbest 



Now velocities are calculated 
using equation 



IfV a+1) <V dmin then V (i+J) = V d „, 
And 

ifV (i+1) >V dmax then V (i+V = V dmax 



New Searching points are calculated using 
the equation X(i+1) = X(i)+V(i+1) 



Check for Constraints, if its not 
violated. Accept it 



Evaluate the fitness values for ne^v searching 
point. If evaluated values of each agent is 
better than precious Pbest, then set to Pbest- if the 
best Pbest Is better than gbesn then set to gbest 




Fig. 2 Detailed flow-chart for applying PSO to RPO 

VI. Results 

The parameters for PSO are taken as per table: 

The anticipated algorithm is run for objective function 
i.e. minimization of real power losses for different cases. 
As main objective is reactive power optimization so we 
deviate generator bus voltages (from 1.0 to 1.1 in case 1 & 
case 2, from 0.9 to 1.0 in case 3 & case 4) and shunt 
capacitor bank setting (from 0.0 to 5.0 in case 1, case 3, 
from 0.0 to 3.0 in case 2, case 4). The consolidated results 
of various cases are shown below in table 

Table: Optimal Parameter Setting for PSO 



Number of iterations 


300 


Cognitive constant, Ci 


2.0 


Social constant, c 2 


2.0 


Max. and Min. inertia 
weights W 


0.4 and 0.95 


Population size 


50 
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Table : Results of IEEE-30 Bus System 



Control 




Case 1 




Case 2 




Case 3 




Case 4 


Parameters 




Constraints 




Constraints 




Constraints 




Constraints 


Gen. Bus 
Voltage (V) 




1 < ^ < 1.1 




1 < ^ < 1.1 




0.9 < ^ < 1.0 




0.9 < \\ < 1.0 


Trans. Tap 
netting yi) 




1.0 <T k < 1.1 




1.0 <T k < 1.1 




1.0 < T k < 1.1 




1.0 <T k < 1.1 


VAR Inst. 

(MVAR) 
\\lc) 




0.0 < Q ci < 5.0 




0.0 < Q ci < 3.0 




0.0 < Q L -i < 5.0 




0.0 < Q ci < 3.0 


Control 




Proposed PSO 
Algorithm 




Proposed PSO 
Algorithm 




Proposed PSO 
Algorithm 




Proposed PSO 
Algorithm 


Parameter 




Min. 


Max. 


Avg. 




Min. 


Max. 


Avg. 




Min. 


Max. 


Avg. 




Min. 


Max. 


Avg. 


v 

Vi 




! 


1.1 


1.053 






1.1 


1.0521 




0.9 


1 


0.953 




0.9 


1 


0.951 


\7 

V2 




! 


1.098 


1.047 






1.1 


1.0455 




0.9 


1 


0.948 




0.901 


1 


0.95 


V5 




1.002 


1.1 


1.049 






1.1 


1.0429 




0.9 


1 


0.951 




0.9 


1 


0.951 


v 




! 


1.1 


1.048 




1.006 


1.1 


1.0576 




0.9 


1 


0.949 




0.9 


1 


0.953 


Vn 




! 


1.1 


1.041 






1.1 


1.052 




0.902 


1 


0.949 




0.901 


1 


0.954 


v 

Vi3 




1.003 


1.099 


1.049 






1.1 


1.0557 




0.9 


0.997 


0.95 




0.9 


0.998 


0.956 


T 

All 






1.1 


1.05 






1.099 


1.0505 




1.004 


1.1 


1.052 




1.001 


1.1 


1.05 


T 

A 12 






1.1 


1.056 




1.001 


1.1 


1.046 




1 


1.098 


1.046 




1 


1.1 


1.058 


T 

A 15 






1.1 


1.048 






1.1 


1.0565 




1.004 


1.099 


1.048 




1 


1.1 


1.05 


T 

A 36 






1.096 


1.047 






1.098 


1.0479 




1 


1.1 


1.043 




1 


1.1 


1.053 






0 


5 


2.728 




0 


2.874 


1.4772 




0.074 


5 


2.838 




0.135 


3 


1.457 






0.052 


4.943 


2.62 




0 


3 


1.5914 




0 


5 


2.041 




0 


3 


1.586 






0 


4.953 


2.223 




0.184 


2.987 


1.5493 




0 


5 


2.822 




0.058 


2.991 


1.771 




































QCl7 




0 


4.706 


2.546 






0 


2.958 


1.3111 




0 


5 


2.649 




0.006 


3 


1.371 


Qc 2 o 




0.068 


5 


2.593 




0.028 


3 


1.5102 




0 


5 


2.553 




0.175 


3 


1.844 


QC 2 1 




0.148 


5 


2.564 






0 


2.915 


1.3708 




0 


4.917 


2.58 




0.026 


3 


1.421 


QC23 




0 


4.997 


2.081 






0 


3 


1.5865 




0.037 


5 


2.579 




0 


3 


1.307 


Qc 24 




0.04 


5 


2.673 




0.014 


3 


1.4088 




0.07 


5 


3.036 




0.088 


2.992 


1.707 


Qc 29 




0.13 


5 


2.436 




0 


2.858 


1.3568 




0.068 


5 


2.878 




0 


3 


1.628 


Power Loss 
(in MW) 




5.173 


7.997 


6.332 




5.289 


7.46 


6.2285 




6.346 


10.47 


7.74 




6.609 


10.18 


7.85 
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VII. Conclusion 
For solving ORPD problem, a new improved integer 
coding Particle Swarm Algorithm is used. The main 
purpose is to minimize the active power loss in the 
network, while satisfying all the power system operation 
constraints. The simulation results show that PSO 
algorithm always leads to a better result. The detailed 
results for different cases have been mentioned in table 
above. As discussed above, the algorithm reaches 
minimum loss under different boundary conditions. After 
studying tables for various cases, the optimum values of 
the control variables are specified in tables in different 
cases of boundary conditions. 
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